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ABSTRACT: 

One and one-half hour isentropic trajectories obtained from an objective 
computer technique for the 308K, 312K, 320K and 332K isentropic surfaces are 
investigated for the May 28, 1967, severe local storm case in the National 
Severe Storm Laboratory (NSSL) network. In conjunction with the trajectories, 
air parcel properties of vertical motion, acceleration, divergence, vorticity, 
static stability and total energy have been studied with respect to development, 
growth, movement and decay of the storm system. Basic assumptions are that 
the flow is frictionless and adiabatic over the entire grid. The results show 
that the observed winds do provide some skill in formulating isentropic trajec- 
tories which demonstrate the presence of storm activity in the specific region. 
The investigation of a mesoscale inertial— gravity wave having a wave length 
of 150 KM, period of 3 hours and phase speed of 14 m sec ^ is also presented, 
with its possible correlation to the onset of the severe storm. 
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ABSTRACT 



One and one-half hour isentropic trajectories obtained from an 
objective computer technique for the 308K, 312K, 320K and 332K isen- 
tropic surfaces are investigated for the May 28, 1967 severe local 
storm case in the National Severe Storms Laboratory (NSSL) network. 
In conjunction with the trajectories, air parcel properties of 
vertical motion, acceleration, divergence, vorticity, static sta- 
bility and total energy have been studied with respect to develop- 
ment, growth, movement and decay of the storm system. Basic 
assumptions are that the flow is frictionless and adiabatic over 
the entire grid. The results show that the observed winds do 
provide some skill in formulating isentropic trajectories which 
demonstrate the presence of storm activity in the specific region. 
The investigation of a mesoscale inertial-gravity wave having a 
wave length of 150 KM, period of 3 hours and phase speed of 14 m 
sec ^ is also presented, with its possible correlation to the 
onset of the severe storm. 
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1. INTRODUCTION 



The severe local storm is an extremely complex meteorological 
phenomenon which can lend itself to interesting and challenging 
research. The approaches one can take for extensive studies are 
almost unlimited as evidenced by the wealth of literature on various 
aspects of the severe storm. Yet there are many questions on the 
formation, movement, growth and decay of these systems which remain 
unanswered. Similar synoptic conditions may seem to exist during 
two separate time periods, and yet storms may be forecasted to 
develop and then not form, or vice-versa. In particular, the extent 
to which an individual storm may evolve is not reasonably understood 
or predictable. 

A. SEVERE STORM RESEARCH 

In 19^+7 a comprehensive study was undertaken in the Thunder- 
storm Project which was perhaps the beginning of extensive scientific 
efforts into the study of severe local storms. From this project 
Byers and Braham (19^9) provided new concepts on the structural fea- 
tures of such storms. The National Severe Storms Project, initiated 
in 1961, was an effort designed to study Great Plains convective 
storms and the environment in which they form. More recently the 
Barbadoes Experiment, which has provided data for study of tropical 
cumulus convection, and Project Storm Fury, for the study of hurri- 
canes, have enabled investigators to carry on detailed research into 
severe storm characteristics. 

Mesoscale networks located in various regions, with instrument- 
ation positioned at specified spacing, have provided for collection 
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of data for certain research objectives. In particular, the National 
Severe Storms Laboratory (NSSL), Norman, Oklahoma, has a systematic 
network of surface and upper air reporting stations and radar facili- 
ties . 

There have been several studies involving kinematic properties 
which have given insights into storm characteristics. Williams (1962) 
investigated small-scale kinematic properties of surface horizontal 
velocity divergence, horizontal stretching deformation, and the 
vertical component of horizontal shearing deformation for isolated 
thunderstorms. He found the intensity and rate of change of these 
properties for small scale systems were 2 to 3 orders of magnitude 
greater than those normally measured on the gross synoptic scale. 
Endlich and Mancuso ( 1968 ) investigated a number of kinematic 
quantities using objective analysis procedures and compared them 
with storm developments. They found vorticity production in low 
levels and horizontal divergence in the upper troposphere to be 
good objective indicators of severe storm activity. Fankhauser 
(1964) used a statistical and observational approach to propose 
that the translation of convective complexes is related in a 
systematic way to the mean wind through the tropospheric layer in 
which the clouds are imbedded, and that this relationship is depen- 
dent on the unique character of the mean wind field. Hammond ( 1967 ) 
investigated an isolated airmass storm, which moved rapidly north- 
eastward across central Oklahoma where its movement was to the left 
of all winds in the middle and upper troposphere. He concluded that 
the same environmental wind field can account for the characteristic 
organization and structure of either right or left moving storms. 
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depending on storm velocity. Studies by Browning and Fujita (1965), 
Newton and Fankhauser (1964) and Browning and Lud lum (1962) have 
shown that large thunderstorms normally travel to the right of the 
winds in the lower and middle level of the troposphere. Matsumoto, 
Ninomiya and Akiyama (1967) studied the relations of cumulus con- 
vections and the mesoscale convergence field over the coastal area 

of Japan and found that well developed cumuli exist only in the area 

-4 -1 

of surface convergence of more than 10 sec 

B. ISENTROPIC ANALYSIS AND TRAJECTORY STUDIES 

The use of isentropic analysis to study the physical processes 
at work in the atmosphere is well known since Rossby (1937) intro- 
duced the technique. The analysis and its uses are well described 
by Namias (1940). Tongues of dry and moist air from specific humidity 
analysis on isentropic surfaces show moisture sources and instability 
regions for thunderstorm activity. 

Calculations of isentropic trajectories for the study of moisture 
sources, vertical motions, velocity accelerations, type of air and 
circulations within systems, vorticity sources, and for investigations 
of transports of pollutants and trace substances has long been appeal- 
ing since non-precipitating motions in the free atmosphere may be 
assumed to be isentropic for limited periods of time. Many investi- 
gators have used graphical calculations from isotach and streamline 
analysis to determine the trajectories (Danielsen, 1959; Reiter 
and Danielsen, 1960; Reiter, 1963). The use of isentropic 
trajectories in thunderstorm studies (Bradley and Feteris , 1966; 
Feteris and Stout, 1966; Carlson and Ludlum, 1968) can provide 
useful information in the study of these mesoscale systems. 
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It is intended in this research to investigate, by the use of 
isentropic trajectories, meteorological properties of air parcels 
and their relationship to severe local storms. Several questions 
will be considered: 

1. Do isentropic trajectories over a short period of time 
describe flow patterns which adequately show the formation, growth, 
movement, and decay of the severe storm? 

2. Do associated parameters computed along isentropic 
trajectories describe motions directly related to the dynamics of 
the severe storm? 

3. Does a combination of any computed parameters correlate 
well with the storm? 

4. Is there an oscillation or wave phenomenon, in particular 
a gravity-inertial wave, that can be detected from the data and/or 
computed parameters which might relate to the local severe storm? 

It is intended to deal with these specific questions in the 
ensuing sections by conducting a case study of a specific severe 
storm system. The isentropic trajectories will be computed from 
an objective computer technique in order that a sufficient number 
of trajectories and isentropic surfaces can be studied for the severe 
storm case. 
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2. OBJECTIVE ISENTROPIC TRAJECTORY COMPUTER MODEL 



Danielsen (1961) introduced an objective method for computing 
isentropic trajectories based on the total energy equation. This 
provided the foundation for the technique for which Mahlman (1968) 
developed a computer method for an objective computation of adia- 
batic three-dimensional trajectories. 



A. BASIC DESIGN 

Utilizing the graphical method by Danielsen (1961), the 
computer program applies an objective method which confines the 
computed trajectory from deviating too far laterally from its true 
position. The concept of the system eminates from the horizontal 
equation of motion in isentropic coordinates, which is 




- f Ik x \V - VM + IF 
0 



(1) 



where \V is the two dimensional horizontal wind vector, f the 

Coriolis parameter, |k the unit vector in the vertical, V the two 

0 

dimensional gradient operator on a surface of constant potential 
temperature 0, M the Montgomery stream function, and |F the friction 
force. 

The basic model utilizes variables of true wind direction, wind 
speed, pressure and Montgomery stream function at each point in a 
four-dimensional (x, y, 0, t) grid. The first step in the program 
is to obtain u and v components of the velocity field, with the 
effects of latitude/longitude type griding included in the calcu- 
lations over the data area. Then, the u and v fields are linearly 
interpolated to give values at twelve time step intervals. 
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The first approximation of the trajectory is using a straight- 
forward kinematic approach. Beginning at a grid point, a parcel is 
displaced for one time step using the initial values of u and v to 
determine the length and direction of the displacement. This dis- 
placement is determined by: 

S = \VAt (2) 

where S is the displacement and \V is the average wind velocity 
along the trajectory over a At interval. Then the time interpolated 
u and v values are spatially interpolated from the four nearest grid 
points to the parcel to obtain new u and v wind components at the 
parcel’s new position. These new values are then used to extrapolate 
the parcel for the next time step. This process is repeated until 
twleve time steps are reached. At each point along the trajectory a 
test is made to determine if the trajectory remains within the original 
grid. If not, then the trajectory is terminated at the point where it 
leaves the field and a new trajectory is begun at the next grid point. 
These isentropic kinematic trajectories can provide a valid result as 
long as the original wind fields have been adequately specified, and 
the time interpolations are satisfactory. 

It is important to note that for inclusion of changes in the 
Montgomery stream function, Mahlman (1968) indicated that the cor- 
rections were sensitive to analysis errors in the M field, rapid 
changes of wind direction perpendicular to the estimated trajectory 
path, and for cases where the wind speed is less than 10 m sec~^. 

Since preliminary studies of the data indicated cases where the wind 
flow was perpendicular to isolines of Montgomery stream function and 
wind speeds less than 10 m sec”"*", the subsequent research was thus 



conducted from the isentropic kinematic trajectories. The incorpo- 
ration of adjustments for changes in Montgomery stream function is 
a separate subroutine in the program so, if only kinematic tra- 
jectories are desired, the subroutine can be left out. The program- 
ming is included in Appendix B. 

B. MODIFICATIONS FOR SEVERE STORM RESEARCH 

An additional field of relative humidity was added in order 
to include a moisture constituent. Routines to compute certain 
meteorological parameters were incorporated so that initial con- 
ditions and variations along the kinematic isentropic trajectories 
could be studied. These routines include calculations of: 




where 5 is the absolute vorticity, D Q the horizontal divergence, 

0 y 

Aq the stretching defor; ation , Bq the shearing deformation, and a 
the static stability. Also the height (z), specific humidity (q), 
saturated specific humidity (q s ) and temperature (T) fields for 
each isentropic surface and data time were determined. Along each 
kinematic trajectory the initial and final values of specific humidity, 
saturated specific humidity, divergence, absolute vorticity, Montgomery 
stream function, kinetic energy, total energy, static stability, and 
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potential vorticity were computed. Total energy, in units of 10 



ergs gm ^ , was defined as: 

TE = C T + gZ + Lq + 

p 2 

where is the specific heat at constant pressure and L the latent 
heat. Potential vorticity was computed as: 

P = °C 0 

where ^ is in units of 10 "^°K sec cm gm \ 

The latitude and longitude positions, u and v components, and 
interpolated height at each of the twelve time steps of the tra- 
jectory were obtained. The average vertical motion (oj) along the 
trajectory was calculated using the procedure: 

- p f p i 
0) = 



( 8 ) 



(9) 



At 



( 10 ) 



-1 



with oj in mb hr and where the i and f refer to initial and final 

pressure over the interval At. The acceleration of the parcel over 

the trajectory was also calculated. Finally, the pressure at which 

the air parcel would become saturated (P ) and the time of saturation 

s 

(t ) were determined by: 



P s = P i[ 1 - 0+ 






SI 

p - p. 

S 1 A J_ 

*s 4t 



T 

P 



( 11 ) 

(12) 



f l 

where the i and f refer to initial and final points on the trajectory 
and P, q and q^ are as previously defined. 

A routine to utilize the CALC0MP plotter with the IBM OS/360 com- 
puter was included in order to graphically display the computed tra- 
jectories on a x, y grid system. The data generated by the program 
was displayed such that subsequent studies could be made. 
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3. DATA AND STORM CHARACTERISTICS 



Data that was considered in the case study are from observations 
in the NSSL mesoscale network on May 28, 1967- Upper air data from 
nine radiosonde stations systematically spaced over the network region 
(Fig. l) provided data at 1700, 1830, 2000, 2130 and 2300 GMT. 

Utilizing computer processing of mesoscale rawinsonde data 
(Kreitzberg and Brockman, 1961 ), fields of Montgomery stream function 
( 10 ^ erg gnT^), wind direction (deg from north), wind speed (m sec -1 ), 
pressure (mb) and relative humidity {%) from the isentropic plane 
section were hand analyzed over the network area. Allowances for 
balloon drift and non-simultaneous releases were not taken into 
account. However, a comparison was made between the original analysis 
and a re-analysis taking into account balloon drift and storm location. 
These results are reported on in later sections. A 9 x 11 latitude- 
longitude grid (Fig. l), with grid separation of .25 degrees, was 
established (southeastern boundary -3^N, 97W). Analyzed data were 
then interpolated to the nearest grid point. Theta surfaces 308 k, 

312K, 320K and 332K were those selected for the analysis levels. 

The 308K theta surface was the lowest level for which data were 
available at most stations over the time periods. The 332K theta 
surface was selected as being representative of the upper level of 
storm influence. This surface was located at a height of about 300 
mb. 

A. BRIEF SYNOPTIC ANALYSIS 

The surface synoptic situation over Oklahoma at 1800 GMT on 
May 28, 1967 indicated a weak high pressure (1018 mb) over the 
Oklahoma Panhaldle, with low pressure areas located over NW Missouri 
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(1010 mb) and in central New Mexico (1005 mb). A weak cold front lay 
between the two low areas, passing through central Oklahoma. The 
front remained nearly stationary during the period under consideration, 
but began to dissipate and become a warm front by 2400 GMT. 

The 1800 and 2400 GMT 500 mb height synoptic analyses showed 
straight southwesterly flow of about 15 kt . The 300 mb height chart 
indicated westerly flow of 20-30 kt at 1200 GMT with winds shifting 
to southwesterly and increasing slightly by 2400 GMT. 

As reported by Fankhauser (1968), the characteristics were only 
marginally conducive to thunderstorm formation, with slightly negative 
stability indexes on morning soundings. The stationary front provided 
sufficient low-level convergence to set off weak convective storms over 
the northwest portion of the network by midday. 

B. STORM CHARACTERISTICS 

The initial radar presentation available for study was at 1817 
GMT. Small echoes were present in the western grid area about 30-50 
nautical miles from station TIK at this time. Larger and more intense 
echoes were present by 2000 GMT. Storm tops to the northwest of 
station TIK exceeded 11 KM in height, and the storm system was oriented 
SW-NE across the grid area. The 2100 GMT radar reflection indicated 
the system to the west of station TIK had decreased in intensity and 
size, with tops reaching to about 9 KM. 

Storm motions used for calculating the relative trajectories 
were obtained from viewing surface radar presentations of similar 
echoes. These velocities were determined as: 
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TIME DIFFERENTIAL COURSE (deg from) SPEED (m sec" 1 ) 



1700-1830 GMT 



0.0 (no echo) 



1830-2000 



268 



1.5 



2000-2130 



2jk 



2.1 



2130-2300 



272 



4.1 



Fankhauser (1968) indicated that echo velocity was 247 (deg 
from) at a speed of 8 m sec Also this course and speed were 

remarkably constant throughout the lifetime of the storm. The 
mean wind direction through the atmospheric layer for the 1700- 
2300 GMT period was 247 (deg from) and had a mean wind speed of 
8 m sec . 
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k. DISCUSSION OF RESULTS 



A. SEVERE STORM RELATIONSHIPS TO BASIC DATA ANALYSES 

There were five basic fields which were analyzed on the isen- 
tropic surfaces in the mesoscale network. No objective smoothing was 
performed on these fields. Part of the research was to investigate 
these basic fields and to determine if they provide an indication of 
storm development, growth, movement and decay. 

The mesoscale analysis over the NSSL network disclosed that a 
region of lower Montgomery stream function values was maintained 
near the middle of the grid area and to the southwest of the storm 
region. This prevailing pattern was consistent on each theta surface. 
Values were generally larger at the initial data time of 1700 GMT and 
then decreased in value over the entire grid as time progressed. How- 
ever, the predominant contour pattern remained the same (values on the 
308K theta surface at 1700 GMT ranged from 310. Uo joules gm to 
310.62; at 2130 GMT values ranged from 310.17 to 310.45). The hori- 
zontal gradient (V M) was stronger to the SE of the storm area. A 

y 

stronger gradient existed at 2300 GMT on the theta surfaces than for 
previous times. It was noted that the thunderstorm had almost dissi- 
pated by this time. A dome of colder temperatures was associated 
with the low stream function values. This low M value and strong 
gradient region could be associated with the stationary front which 
existed over central Oklahoma and would indicate the thunderstorm was 
occurring northwest of the front. Figure 2a is an example of the 
Montgomery stream function pattern. 
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Figure 2a, 320K Theta surface# 2000 GMT Montgomery stream 
function analysis (add 320,00 to each value ), Heavy solid 
line represents region of storm activity. Values in joules 
gm-'. 
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conducted from the isentropic kinematic trajectories. The incorpo- 
ration of adjustments for changes in Montgomery stream function is 
a separate subroutine in the program so, if only kinematic tra- 
jectories are desired, the subroutine can be left out. The program- 
ming is included in Appendix B. 



B. MODIFICATIONS FOR SEVERE STORM RESEARCH 



An additional field of relative humidity was added in order 
to include a moisture constituent. Routines to compute certain 
meteorological parameters were incorporated so that initial con- 
ditions and variations along the kinematic isentropic trajectories 
could be studied. These routines include calculations of: 



D e = 



A e = 



B e = 



o r 



_ ftv _ 9u\ 

v* 9y/ , 

— 4 

x / 3y / 1 



+ f 




(3) 

(U) 

(5) 

( 6 ) 
( 7 ) 



where c is the absolute vorticity, D fl the horizontal divergence, 

Aq the stretching deformation, B^ the shedding deformation, and a 

the static stability. Also the height (z), specific humidity (q), 

\ 

saturated specific humidity (q s ) and temperature (T) fields for 

each isentropic surface and data time were determined. Along each 

\ 

kinematic trajectory the initial and final values of specific humidity, 
saturated specific humidity, divergence, absolute vorticity, Montgomery 
stream function, kinetic energy, total energy, static stability, and 
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Pressure analyses on each theta surface exhibited sporadic 



patterns. There was no real evidence of vertical consistency with 
regions of high or low pressure, and no time dependency as to move- 
ment or preservation of any defined center. There were large pres- 
sure ranges in the grid area on the 308K and 312K theta surfaces at 
2130 and 2300 GMT ( 7 U 8-865 mb for 308K and 649-753 mb for 312K). 

These pressure ranges were a result of the decrease in height of the 
308K and 312K surfaces in the SW corner of the grid at 2130 GMT and 
extending into the center of the grid by 2300 GMT. Pressure analyses 
on the 320K theta surface shows that a strong gradient existed to the 
SE of the storm. An example of the pressure pattern is given in 
Figure 2b. 

An analysis of relative humidity revealed a region in the S-SE 
area of the network which was a minimum value area on all theta sur- 
faces at all times (Fig. 2c). Since the storm activity occurred 
between upper air network stations, there was no clear indication of 
a region of high relative humidity. A moist or dry tongue was not 
evident from the relative humidity analyses. 

There was a region of higher wind speed in the NW corner of the 
grid area and lower wind speeds in the SE corner (Fig. 2d) on all 
theta surfaces. This wind maximum propagated toward the east with 
time and increased in intensity until 2300 GMT when the speed began 
to decrease. 

Wind direction was generally anticyclonic in the SE area and 
cyclonic in the NW region with a shear line along the thunderstorm 
activity (Fig. 2e). Time sections at stations LTS, FSI, CHK and TIK, 
along with the theta surface analyses of wind direction, indicated 
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Figure 2b, 320K Theta surface* 2000 GI-TT preaaure (mb) analysis^ 

Haary aolid line represent* region of atom activity. 
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Figure 2c, 320K Theta surface# 2000 GMT relative humidity 
(%) analyst*. Heavy solid line represent* region of *tonn 
activity. 
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Figure 2d, 320K Theta surface, 2000 GMT wind speed (m sec"' ) 
analysis, Heavy solid line represents region of storm activity. 
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Figure 2m, 320K Thata surfaca, 2000 GMT wind direction (dags 

from north) analysis. Heavy solid line represents region of 
storm activity. 
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the wind veered with height and with time on the 308K, 312K and 320K 
theta surfaces. The time sections indicated the wind began to back 
with height above the 320K surface. This backing was evident by 1830 
GMT at station LTS , by 2000 GMT at FSI, by 2130 GMT at CHK and by 
2130 GMT at TIK. Wind speed at these stations generally increased 
with time and with height. 

Data from stations LTS and FSI, which are toward the SW of the 
storm area, indicated the low level wind speed increased until about 
2015 GMT and then began to decrease. Stations CHK and TIK, located 
in the SE edge of the storm region, show the wind speed continues to 
increase over the whole time period (1700-2300 GMT). Both the veering 
with height and speed increase have been correlated with severe storm 
activity (Miller, 1967 ) . 

B. TRAJECTORY RELATIONSHIPS TO THE SEVERE STORM 

The kinematic isentropic trajectories were computed for one and 
one-half hour time periods at each grid point and on each theta surface. 
For reference to the trajectories that are discussed, see Appendix A. 

The 1700-1830 GMT trajectories clearly show veering with height 
and time in the wind field. The 308K and 312K trajectories indicate 
anticyclonic flow over the whole grid, whereas the 320K and 332K 
trajectories have cyclonic flow in the W-NW corner and anticyclonic 
flow in the S-SE area. This is the flow pattern just prior to the 
storm formation, since the thunderstorm did not appear on radar until 
about 1817 GMT. 

The 1830-2000 GMT trajectories imply changes in the wind field 
that are significant. On the 308K surface, the trajectories are 
beginning to have cyclonic flow in the W-NW, with a shear line almost 
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along the line of thunderstorm activity. The 312K trajectories specify 
anticyclonic flow over the whole area. The 320K surface flow pattern 
exhibits cyclonic flow along the south edge of the storm area, anti- 
cyclonic flow to the north, and straight flow along the center line 
of the storm activity. The 332K trajectory pattern denotes cyclonic 
flow to the SW and along the southern edge of storm activity, with 
anticyclonic flow in the NW and SE grid area. The storm appears to be 
oriented along the 320K trajectory flow pattern over this 1830-2000 GMT 
time frame. 

During the time interval 2000-2130 GMT, the storm was in its most 
intense stages. It appears to be more oriented along the 308K flow 
pattern, which is from the south and cyclonic over most of the grid. 

The area to the west of the storm is strongly cyclonic. The 312K 
trajectories indicate straight SW flow. The 320K trajectories give 
evidence of a cyclonic area just to the south of the more intense 
regions of the storm. The 332K trajectories display anticyclonic flow 
over the entire grid. There is a veering of 70 degrees or more between 
the 308K and 320K trajectories over this period and pronounced backing 
in the NW region between the 320K and 332K surfaces. 

Radar presentations indicate the dissipation of the storm over the 
2130-2300 GMT period. The 312K trajectories are from the SW and show 
slight anticyclonic turning. The 320K and 332K trajectories indicate 
a slight cyclonic turning in the western portion of the grid and anti- 
cyclonic turning in the eastern portion. Strong backing of degrees 
or more with height is evident between the 320K and 312K trajectory 
patterns . 

Selected composite trajectories were obtained from each time interval 
in order to investigate air parcel flow over the whole time period . 
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These trajectories and associated parameters are discussed in the 
following section. 

C. RELATIONSHIP OF ASSOCIATED PARAMETERS TO THE SEVERE STORM 

One of the most important parameters to he computed from the 
kinematic isentropic trajectory program was the average vertical motion 
distribution for each theta surface and each 1.5 hour time interval. 

When used alone, these average vertical motion patterns are not suffi- 
cient to describe the storm activity. However, they do provide very 
definite clues to the storm dynamics when used in conjunction with 
other properties. 

During the time interval 1700-1830 GMT, the vertical distribution of 
vertical motions over the NSSL network indicated lower level air parcel 
sinking in the western half of the grid of about 8 mb hr ^ and rising 
motions of 10-20 mb hr ^ in the eastern half of the grid. Upper level 
patterns were reversed, with rising motions in the western half and 
sinking motions in the eastern half. In addition, there was very strong 
rising motion (20-30 mb hr in the western half of the grid at a height 
of about 10 KM. It is possible that the high level vertical motion pat- 
tern is the initial indication of an upper level formation of supporting 
vertical motions for the storm. 

The 1830-2000 GMT vertical motion patterns (Fig. 3a-3d) clearly 
delineate and support the storm system. At each level there is strong 
rising motion as air parcels come into the storm area, with sinking 
motions over the remaining grid area. The 320K isentropic surface 
vertical motion pattern indicates the region of maximum rising motion. 

The storm is also oriented parallel to the 320K trajectory flow pattern 
over this period. It would appear the 320K (i.e. approximately 500 mb) 
flow is responsible for maintaining the storm at this time. 
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Over the 2000-2130 GMT period the 308K and 312K isentropic sur- 
faces vertical motion patterns present sinking motion over most of 
the grid, with larger values (25-30 mb hr in the northern half of 
the grid. The 320K and 332K theta surfaces vertical motions still 
show rising in front and into the storm system, but the patterns are 
not as well defined. The 2000-2130 GMT convergence fields on the 320K 

isentropic surface are well correlated with the storm system and show 

-4 -1 

strong convergence of 1.0 to 1.5 x 10 sec in the system (Figs. 4a, 
4b). 

Darkow (1968) has examined the total energy environment of severe 

storms on the synoptic scale of motion and has determined a "Total 

Energy Index”. This index is computed as the difference between the 

"total energy” (c T + gz + Lq) at 500 mb and the total energy at 850 
P 

mb. Index values more negative than -2.0 cal gm indicate severe storms 
and associated tornado activity are highly probable. 

The NSSL isentropic plane section computer output for the nine 
radiosonde stations gives the static energy. This data was investi- 
gated in terms of a "Total Energy Index”. Values at 1700 GMT 
indicated a region of -2.0 in the area where the thunderstorms 
began to occur. Values in the same region were -4.0 or larger 
by 2000 GMT (Fig. 5). Values over the remaining grid area were 
slightly negative or positive. 

Total energy, defined by equation (8), was investigated for 
storm evolvement. The 312K and 320K isentropic surfaces for the 2000 
and 2130 GMT analyses show a region of maximum total energy in the storm 
vicinity, with decreasing values toward the NW and SE. This region is 
seen in Figure 6. As shown by Figure 7, the contour pattern of total 
energy closely resembles the specific humidity distribution, which 
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Figure 4a, 320K. Thata aurfaca, 2000 GMT convergence (10~ % aac“') 
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Figure 5, 2000 GMT Total energy index. 
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indicates the Lq term dominates the pattern. Thus the relative humidity 
analysis is a critical factor in determining the total energy. The 
relative humidity is not the most reliable information from a radio- 
sonde sounding, particularly at higher altitudes. Since the Lq term 
seems to dominate the total energy contours, unless the relative 
humidity field is adequately specified, the total energy computations 
will not provide good information with which to study the storm. There 
axe also inaccuracies in the Montgomery stream function field which may 
amplify the error in total energy computations. 

Fankhauser ( 1968 ) has demonstrated the wind field relationship for 
this same May 28, 1967 storm system. Using chaff trajectories and 
aircraft measurements, the streamline and isotach pattern of the com- 
posite relative wind field with respect to the storm echo are presented 
(Fig. 8 ). The 320 K isentropic surface almost corresponds with the 500 
mb isobaric surface, thus a good comparison can be made. In Figure 9 
the isentropic trajectories show that flow is being deviated around the 
storm, convergence into the storm and divergence behind the storm. The 
acceleration field (Fig. 10) also shows maxima to the right and left of 
the storm and minima in front of and behind the storm. Fankhauser has 
demonstrated the flow for an isolated storm echo and the fields are 
instantaneous. The isentropic trajectories show average conditions 
for the larger storm system, and yet the results are strikingly similar. 
The comparison times are also different, which may indicate this wind 
distribution is present throughout the storm cycle for both isolated 
cells and the entire storm complex. Thus the 1.5 hour kinematic isen- 
tropic trajectories correlate very well with the chaff trajectories 
that Fankhauser has obtained. 
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Figure 8 , 2145 GMT simplified streamline and isotach pattern 

relating the composite relative wind field with respect to 
storm echo centered on cross, (after Fankhauser# 1968 # fig, 8 ), 
Wind speeds in m sec"'. 




Figure 9, 320K Theta surface# 1630-2000 GMT kinematic 

isentropic trajacteriaa relative to the atom. Dashed 
limes represent changes in absolute vorticity (lO^sec*) 
computed along the trajectory. Heavy solid lines represent 
changes in divergence (lO^sec* 1 ) computed along the 
trajectory. Storm position as of 1915 GMT, 
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Figure 10, 320K Theta surface, 1830-2000 GMT acceleration* 

(a sec*' hr" 1 ). Storm position as of 1915 GMT, 



D. INERTIAL-GRAVITY WAVE ASSOCIATION TO THE SEVERE STORM 



A close look at various analyses seems to imply an oscillation 
present in the network area. Thus an investigation into the possibility 
of a wave like phenomenon present over the data area was considered as 
a possible source of the storm activity and storm movement. 

An oscillatory tendency can be observed on the 320K theta surface 
composite trajectories (Fig. ll). Measurements in the center of the 
grid region show a wave length of about 150 KM and period of about 3 
hours. Radar presentations of the most intense storm activity also 
indicate a wave length of about 150 KM (Fig. 12). 

Surface pressure tendencies over the network area were studied to 
determine if the oscillation could be detected at the surface. Micro- 
barograph traces were available at 35 stations over the data area. The 
pressure tendencies indicated falling pressures of about 0.03 inches of 
Hg per 1.5 hours over the grid during the entire life cycle of the storm. 
Figure 13 illustrates this tendency over the area. A noticeable fluctu- 
ation was not apparent at most stations. Even the presence of the storm 
was difficult to detect (see stations 3A, 3C, UB and Uc on Fig. 13). 

Synoptic pressure tendency patterns were obtained for 1800 and 2^00 
GMT, and a mean synoptic pressure tendency over the grid area was cal- 
culated. This value was then subtracted from the network stations 
pressure tendency at each time frame in order to calculate a "meso- 
scale pressure tendency 11 . The results, as seen in Figure l4, indicate 
both rise and fall areas. The rise and fall areas show a half wave 
length of about 100 KM along the mean low level wind field, and a 150 
KM wave length between two negative centers along a line of mean wind 
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Time sections at selected stations along the mean wind direction 
(see Fig. 15) were made in order to study the time variations of 
vertical motion, static stability, relative humidity, wind speed and 
wind direction. The time sections at stations LTS, FSI, CHK and TIK 
(Figs. l6a-l6d) indicate an apparent periodicity in the vertical motion 
and static stability patterns. The vertical motions are the average 
omega obtained by the trajectory method, and the static stability is 
the value centered between two adjacent theta surfaces. Measurements 
on these time sections, from the vertical motion patterns, give a wave 

length of about 150 KM, wave period of 3 hours and wave speed of lU m 

-1 

sec 

It can be seen in Figure l6a that a region of vertical motion of 
+30 mb hr ^ is occurring at station LTS at 1915 GMT, between the 308K 
and 312K theta surfaces. With the measured wave speed of Ik m sec ^ 

(50 KM hr 1), a propagating wave without change in amplitude, constant 
wave length and moving in the direction of the mean wind would produce 
a region of about +30 mb hr 1 at station FSI about 20^5 GMT and station 
CHK about 21^5 GMT. This condition is verified at these stations for 
the specified times (Figs. l6b & l6c ) . 

With the assumption of a wave-like phenomenon present, a study was 
conducted to determine if such a wave phenomenon and other supporting 
dynamic and thermodynamic forces could be correlated with the onset of 
the storm. The earliest radar presentation at 1817 GMT indicated small 
cumulus activity in the western portion of the grid area. Time sections 
at stations LTS, FSI, CHK and TIK (Figs. 17a-17b) show that at 1700 GMT 
there is a pronounced shift in wind direction between the stations on 
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Figure 15, Station cross-section distribution for 
time sections. 
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Figure 16a, Tima section at station LTS, Heavy solid lines 
represent vertical motions (mb hr"*). Dashed lines represent 
static stability (K mb"'). Thin solid lines represent theta 
surface positions. 
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Figure 16b, Tin* section at statioa F3I, For •xpl*a«tioa 
fig, 16*, 
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Figure 16c, Time section at station CHK, For explanation see 
fig, 16a, 
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Figure 17a, Time section at station LTS, Thin solid lines ars 
theta surface positions. Heavy solid lines are relative humidity 
(%), Dashed lines represent wind speed (m sec - *). Short solid 
lines represent wind direction (360 degs at top of page). 
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Figure 17c, Tim* section at station CHK, For explanation see 
fig, 17a, 
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Figure 17d, Time section at station TIK, For explanation see 
ti<3, 17a, 
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the 332K theta surface. This wind shift could be further evidence of 
the passage of a wave through the area around 1700 GMT. 

The relationship between the horizontal curvature of the trajectory 
(K^) and the streamline (K^) can be related by Blaton's equation: 

(13) 



K - K + 1_ Iff. 

^ K s V 3t 
s 



3 f f 

where V is in the wind speed and — — is the local rate of turning of the 
S (j *c 

3 f f 

horizontal wind direction. At station FSI (Fig. 17b) — — is positive, 

3 1 

thus should be greater than K , which is what is observed (Fig. 18 ). 

3 •pf' 

At station CHK (Fig. 17c) — is negative, thus should be less than 
K g , which is also the case (Fig. l8). If the streamlines are wave- 
shaped and if the wave moves with a speed greater than the mean wind 
speed, the curvature of the streamlines and trajectories may be of 
opposite sign and the trajectories will have a curvature greater than 
the streamlines . This approach implies that a wave having a phase 
speed greater than the mean wind speed is indeed present in the region. 

The 1700-1830 GMT vertical motion, 1700 GMT static stability and 
divergence patterns were superimposed upon the 1700 GMT streamline flow 
on the 332K theta surface (Fig. 18). A parcel of air, with respect to 
a wave that is traveling faster than the mean wind speed, is actually 
moving toward the west. Thus the air parcels are experiencing con- 
vergence, sinking motions and entering a stability minimum area in the 
western portion of the grid. 

From the conservation of potential vorticity on an isentropic 
surface, which is given by: 



dt U Q o) = 0 



(lU) 
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where 5 and c are defined as in equations (3) and (7)> an air parcel 
entering a trough region would normally experience an increase in 
vorticity, thus requiring a decrease in static stability. From the 
continuity equation in isentropic coordinates in the form: 

at ' °V V = 0 (15) 

it is seen that, where o is positive, a parcel experiencing convergence 
must have a decrease in static stability. 

There thus appears good evidence of the presence of a wave 
phenomenon in the region in the early stages of storm activity. There 
also is a good indication that the wave has a phase speed greater than 
the mean wind speed. The 1700 GMT model on the 332K isentropic surface 
gives evidence that the parameters can only be consistent for a wave 
speed greater than the mean wind speed. With this dynamically consistent 
model there has been established a region conducive to convective 
activity for the 1700-1830 GMT period which could initiate the severe 
storm activity. 

For adiabatic, frictionless flow along a streamline, the equation 
of motion in isentropic coordinates reduces to the form: 



dV 

s 

dt 



as 



(16) 



By assuming: 



at as 

where C is the phase speed, equations l6 and 17 can be reduced to the 
form: 

aM 

s_ 

as 



(17) 



(18) 
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where V g is the mean wind speed. Now by assuming a phase speed of 
twice the mean wind speed, a relationship for adjusting the mass 
field to the wind field can be obtained in the form: 

dV 



cfri T7 o 

3s " V s 3s~ 



( 19 ) 



such that AM a AV g along the streamline for the same As. Figure 19 
illustrates the adjustment that is necessary in the Montgomery stream 
function field. This adjusted field is dynamically consistent with 
the model that has been established on the 332K theta surface at 1700 
GMT, Air parcels, moving with respect to the wave, which are experi- 
encing an acceleration are also experiencing a decrease in Montgomery 
stream function. Thus knowledge of the wave in the region can be 
helpful in adjusting the mass and wind fields such that they are 
dynamica lly consistent . 
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5. PROBLEMS ENCOUNTERED 



One of the most demanding and important parts of the research are 
the initial analyzed fields. The problem has required the analysis of 
100 data fields. Proper analysis involves both spatial and time con- 
siderations, as well as dynamical and thermodynamical relationships 
between variables. The horizontal gradients are very important as 
interpolations are made between network stations. Subjective smoothing 
can destroy the very detail that may well be important , but at the 
same time under smoothing may produce anomalies that invalidate sub- 
sequent studies. Station bias must be carefully considered. 

On the synoptic scale, the normal concern of changes in a certain 
parameter that must be detected over a specified distance or time will 
usually be expressed in, say, whole numbers. However, for mesoscale 
work one must detect changes that amount to tenths or hundredths for a 
specified distance or time. The instrumentation for data collecting is 
the same for both scales of motion, thus the sensitivity of the instru- 
ments and the magnitude of instrument error and measurements is very 
critical . 

Engelmann and Davis (1968) have found that the standard observational 

•7 

error in radiosonde data will give a computational error of 0.0U x 10 
2 -2 

cm sec in the Montgomery stream function. These errors result from 

errors in the instruments and receiving equipment, baseline check and 

in the mean temperature. The M fields were analyzed for each 0.05 x 
7 2-2 

10 cm sec interval, which is about the value of the computational 
error. Thus the M values and analyses fields contain inaccuracies 
which effect calculations using these fields. It was seen in section 
hD that the Montgomery stream function fields required rather large 
adjustments in order to be dynamically consistent with the winds. 
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Figures 20a-20e indicate the changes in the analyses for the 
320K isentropic surface at time 2000 GMT when balloon displacement, 
location of the storm at analysis time, and a conscientious effort 
toward spatial and time continuity were made. The 320K surface was 
selected as a comparison for the re-analysis since its height fluctu- 
ations were minimal over the entire grid and between each data time. 
Isolines of Montgomery stream function were oriented along streamlines 
and rapid changes perpendicular to the flow were minimized. Pressure 
isolines were also oriented along the flow, as were relative humidity 
values. Also larger values of relative humidity were analyzed for, 
in the vicinity of the storm. Wind direction and speed considerations 
were carefully matched, particularly with respect to time continuity. 

The effects of the re-analyses on the trajectories and values of 
the associated parameters were investigated to determine if any sig- 
nificant changes resulted. Figure 21 shows the comparison between 
the initial analysis (IA) and the re-analysis (RA) for selected tra- 
jectories and their associated parameters. Vertical motion (VM) 

(mb hr acceleration (A) (m sec” 1 ), vorticity change (V CHG) (sec” 1 ), 
divergence change (D CHG) (sec 1 ), total energy change (TE CHG) 

(10 T erg gm 1 ), static stability change (SS CHG) (°K mb” 1 ) and potential 
vorticity change (PV CHG) (10 ^ °K sec cm gm ^ are given for each of 
the trajectories. Vertical motions and accelerations are smaller for 
the re-analyzed case and there are two cases which have sign changes. 

The re-analyzed trajectories have more cyclonic curvature in the vicinity 
of the storm, and trajectory 1 has more anticyclonic curvature. Other 
parameters show little difference. 
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Figur* 20a, 320X Thata aurfao*, 2000 GMT Montgomery stream 

fuactio* 9am fig, 2a for axpliaatioa. 
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Figure 20b, 320K Theta surface# 2000 GMT pressure (rob) 

re-analysis. 
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Figttr* 20c, 320K That* iurftc«< 2000 GKT ralatlT* humidity 

WU ra-analysia. 




Flgur* 20d, 320K Th*t* rurf*c«» 2000 SMT triad apaad (a mc~*) 
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Figure 20», 31 OK That* mrfaoe# 2000K vied direct iom ( dmga 

tram aerth) re-«a*lyal«. 
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Using equation 16 and the 2000-2130 GMT accelerations that were 

computed from the kinematic isentropic trajectories, the expected 

changes in the Montgomery stream function values along a selected 

AS of 50 KM were computed at various areas over the grid on the 320K 

theta surface (Fig. 22). The solid black lines are the analyzed M 

field based on the data at the nine NSSL upper air stations at 2000 

7 2 -2 

GMT. The values are given in units of 10 cm sec . The dashed 
black lines represent the 2000-2130 GMT accelerations (m sec*' 1 ' hr”' 1 ') 
and the short green lines are the position of AS over the area. In 
Table 1 a comparison of the predicted value and measured value for 
the M values are given. In many cases there is almost an order of 
magnitude difference, with the predicted values having a much lower 
absolute value. Using these predicted changes and approximate values 
at each station, an adjustment to the M field was made (Fig. 22 - red 
dashed lines). This new M field is now dynamically consistent with 
the winds. It would thus seem that the Montgomery stream function 
and hence the temperature and height fields are in error. It is not 
conclusive however, but along with the adjustments in the M field 
discussed in Section l+D there appears to be unacceptable errors in 
the initial data fields of Montgomery stream function. From this 
point of view, therefore it appears the wind data is more accurate. 
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Figure 22, 320K Theta surface# 2000 GKT Moatgooery stream fuactioa adjustmeat 
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6. COMPARISON WITH THEORY 



In the basic assumption of frictionless, adiabatic flow along 
an isentropic surface, on which the study was based, there are 
three parameters which should be conserved for an air parcel moving 
in the flow. These parameters are potential temperature (©) (°K), 
specific humidity (q) (gm gm and potential vorticity ( ?) (10 10 °K 
sec cm gm “*■). The initial input fields into the computer program were 
analyzed on constant theta surfaces, thus the program constrained theta 
to be conserved for each trajectory. However, specific humidity and 
potential vorticity were allowed to vary. 

Selected composite trajectories were studied in order to determine 
both short (1.5 hour) and long (6 hour) term variations. In Figure 22 
and Table 2 the deviations are shown for the 320K theta surface. 
Vertical motions (U) (mb hr ^), accelerations (A) (m sec ^ hr ^), 
static stability (<T) (°K mb ^) and total energy (TE) (10^ erg gm 1 ) 
were also studied in conjunction with the conserved properties to 
determine if any correlations exist. Trajectories labeled 1A-1F are 
for 1700-1830 GMT; 2A-2F for 1830-2000 GMT; 3A-3F for 2000-2130 GMT; 
and UA-l+F for 2130-2300 GMT. Both initial (i) and final (f) values 
for each of the properties are given. 

With the exception of 1A-1F trajectories, the specific humidity 
was somewhat conserved for each trajectory. Changes are random and 
do not appear to be associated with vertical motions or accelerations. 
Trajectories 2E and 3F indicate their deviations may be from the storm 
effects, since their initial point is outside the storm and final point 
inside the storm region. 
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7- CONCLUSIONS 



The research has demonstrated the critical nature of mesoscale data 
both in instrumentation and analysis errors. It was found that the wind 
data were perhaps the most reliable for use in the study. Subsequently 
the 1.5 hour kinematic isentropic trajectory relationship to the severe 
local storm warranted considerable exploration. It was further noted 
that the 320K theta surface (approximately 500 mb) showed the tra- 
jectory and air parcel properties which were more closely associated 
with the life cycle of the storm. 

The trajectories indicated that prior to storm development the wind 
veered with height at all levels, flow was anticyclonic with the ex- 
ception of cyclonic flow in the NW corner in upper levels, and a wind 
maximum region existed in the NW corner of the network. With the 
presence of an inertial-gravity wave at the 332K theta surface at 1700 
GMT, the dynamically consistent vertical motions, static stability, 
convergence, trajectories, streamlines and adjusted Montgomery stream 
function fields support conditions which could initiate convective 
activity. It was also noted that the storm is to the NE of the 
Wichita mountains, and the existance of a lee wave supporting the 
inertial-gravity wave is worthy of additional studies. 

The study of various air parcel properties on the isentropic 
surfaces has provided insights into the storm dynamics. The 1830- 
2000 GMT vertical motions indicated supporting rising motions for 
storm growth. The 2000-2130 GMT, 320K theta surface convergence 
patterns help to illustrate the intensity of the storm during this 
period. The lack of supporting vertical motions, and the anti- 
cyclonic flow in each level tended to provide the response necessary 
for decay of the storm over the 2130-2300 GMT period. 
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The orientation of the storm along the upper level trajectories 
in the early stages and along the low level trajectories in the latter 
stages may he indicative of storm movement. When the storm orientation 
is considered along with the inertial-gravity wave, there exists the 
possibility the storm formed at high levels and propagated downward 
into the lower levels. It is concluded that the kinematic isentropic 
trajectories, over short periods of time, do provide an effective 
method for studying local severe storm activity in the mesoscale range. 

8. SUGGESTIONS FOR FUTURE DEVELOPMENTS AND IMPROVEMENTS 

Since the initial data fields are critical to the success of the 
computer model, it would be of benefit to investigate the application 
of error analysis techniques in order to improve the initial data. 

Further studies are necessary when corrections for Montgomery stream 
function changes can be made, in order to improve the final trajectories. 
Incorporation of diabatic effects in regions where air parcels are be- 
coming saturated prior to the final trajectory point can be included in 
the model. Further studies of different storm systems is encouraged in 
order to test the short period kinematic isentropic trajectories, and 
to determine if they can be utilized as a forecasting aid. Since there 
existed deviations in conservation of specific humidity and potential 
vorticity along the trajectories, additional research will be necessary 
to investigate the conservation of these properties. With the sup- 
porting evidence of an inertial-gravity wave in the region, a theoreti- 
cal approach is recommended in order to study the phenomenon in the 
range of wave length and wave speed encountered in this case study. 
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APPENDIX A - ISENTROPIC KINEMATIC TRAJECTORIES 



This appendix includes the kinematic trajectories which have been 
computed for the 308K, 312K, 320K and 332K isentropic surfaces for 
the time periods 1700-1830, 1830-2000, 2000-2130 and 2130-2300 GMT. 
They were produced by the objective computer model for the May 28, 
1967 severe storm case. The grid is an x, y grid, with equal grid 
spacing. 
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APPENDIX B - ISENTROPIC TRAJECTORY COMPUTER PROGRAM 



This appendix includes the objective computer model. The program 
is written in FORTRAN language. Each routine contains comment cards 
which explain the routine and each step. 
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